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ABSTRACT 


The  sulfate  diffusion  coefficient  is  a variable  in  theoretical 
models  used  for  predicting  the  service  life  of  concrete.  Analysis 
of  sulfate  models  predicts  that  if  the  process  is  either 
diffusion  or  reaction  controlled,  a 10-fold  increase  in  the 
diffusion  coefficient  results  in  a 10-fold  increase  in  the 
predicted  service  life.  Therefore,  reliable  estimates  of  sulfate 
diffusion  coefficients  are  needed  to  make  reasonable  service  life 
predictions . 

Methods  for  determining  the  sulfate  diffusion  coefficient  of 
sulfates  in  cement-based  materials  are  explored  in  this  paper, 
including  experimental  methods  and  mathematical  modeling.  Also, 
an  approach  for  estimating  the  diffusion  coefficient  of  sulfate 
ions  is  presented,  which  is  based  on  comparing  sulfate  diffusion 
coefficients  with  chloride  diffusion  coefficients. 


Key  words:  building  technology;  concrete;  diffusion  coefficient; 
radioactive  waste;  service  life;  sulfate  attack 
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1.  INTRODUCTION 


Sulfate  attack  is  one  of  the  possible  major  degradation  processes 
for  concrete  vaults  buried  underground  for  the  disposal  of  low- 
level  radioactive  waste  (LLW)  [1].  The  diffusion  coefficient  for 
sulfates  in  the  pore  water  of  concrete  is  a variable  in  most 
models  for  predicting  the  service  life  of  concrete  controlled  by 
sulfate  attack  [2,3].  These  service  life  predictions  are 
sensitive  to  the  values  of  the  sulfate  diffusion  coefficient  [3]. 
Because  of  the  lack  of  data  on  sulfate  diffusion  coefficients, 
the  diffusion  coefficient  for  chloride  ions  in  cement  paste  has 
been  used  as  the  value  of  the  sulfate  diffusion  coefficient  in 
making  service  life  predictions  [3].  Therefore,  information  is 
needed  on  values  of  sulfate  diffusion  coefficient  to  make 
reasonable  service  life  predictions  when  sulfate  is  the 
controlling  degradation  process.  Methods  for  determining  the 
sulfate  diffusion  coefficient  in  cement-based  materials  are 
explored  in  this  paper,  including  experimental  methods  and 
mathematical  modeling.  Also,  a method  for  estimating  the 
diffusion  coefficient  of  sulfate  ions  is  presented,  which  is 
based  on  comparing  sulfate  diffusion  coefficients  with  chloride 
diffusion  coefficients. 

We  are  using  Tntrinsic  diffusion  coefficient,  Di,  to  denote  the 
diffusion  coefficient  in  a porous  medium  without  absorption  or 
reaction  occurring,  while  effective  diffusion  coefficient,  Dg, 
denotes  the  diffusion  coefficient  in  a porous  medium  in  which 
either  or  both  absorption  and  reaction  occurs. 


2.  EXPERIMENTAL  METHODS  FOR  OBTAINING  SULFATE  DIFFUSION 
COEFFICIENTS 

2.1  Determination  of  Sulfate  Diffusion  Coefficient  Using 
Diffusion  Cells. 

Using  a conventional  diffusion  cell,  similar  to  that  described  by 
Clifton,  et  al.  [4],  for  measuring  the  diffusion  coefficient  of 
sulfate  ions  through  cement-based  materials  is  significantly  more 
difficult  using  it  for  measuring  chloride  ion  diffusion 
coefficients.  For  example,  the  concentration  (activity)  of 
chloride  ions  which  have  diffused  through  a specimen  can  be 
monitored  with  an  ion  selective  electrode.  However,  an  ion 
selective  electrode  is  not  available  for  sulfate  ions.  The 
concentration  of  sulfate  can  be  measured  by  several  other  methods 
including  precipitating  the  sulfate  as  barium  sulfate  and 
gravimetrically  determining  the  sulfate  content  [5];  by 
potentiometric  titration  using  a lead  salt  [6];  or  by  a 
turbidimetric  method  described  by  Atkinson,  et  al.  [7].  These 
methods  are  significantly  more  time  consuming  than  direct 
measurements  using  a selective  ion  electrode. 
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Calcium  sulfate,  usually  in  the  form  of  gypsum,  is  blended  with 
clinker  during  grinding.  This  calcium  sulfate  is  added  to  reduce 
the  initial  high  rate  of  rapid  hydration  of  tricalcium  aluminate 
and  prevent  flash  set  caused  by  [8].  Thus,  some  of  the  sulfate 
accumulating  in  the  diffusion  cell  would  be  internal  sulfate  and, 
thus,  would  increase  the  measured  or  effective  diffusion 
coefficient.  Another  complication  in  measuring  sulfate 
diffusivity  using  a diffusion  cell  is  the  possible  degradation  of 
specimens  caused  by  sulfate  reactions  with  constituents  of  the 
cement.  For  example,  diffusing  sulfate  ions  could  react  with 
tricalcium  aluminate  hydration  products  to  form  ettringite. 
Ettringite  occupies  more  volume  than  the  solid  reactants, 
therefore,  its  formation  could  result  in  cracking  of  specimens. 

2.2  Determination  of  Sulfate  Diffusion  Coefficient  by 
Radioactive  Tracers 

A radioactive  traced  method  was  developed  by  Spinks,  et  al.  [9] 
for  studying  the  diffusion  process  of  anions  and  cations  in 
hardened  portland  cement  paste.  They  labeled  salts  with  Ca, 

3^S,  and  ^^Na,  and  exposed  cement  and  mortar  specimens  to 

aqueous  solutions  of  the  labeled  salts.  The  depth  of  diffusion  of 
salts  was  measured  by  cutting  slices  from  the  specimens  and 
measuring  the  radioactivity  of  the  slices  with  a Geiger  counter. 
Effective  diffusion  coefficients  were  calculated  using  Pick's 
first  law  of  diffusion.  For  a 1:2  mortar  (one  part  by  mass  of 
cement  to  2 parts  by  mass  of  Ottawa  sand)  , De  for  sulfate  was 
6.1x10"^^  m^/s,  and  for  iodine  it  was  130x10*^^  m^/s,  giving  a ratio 
for  the  diffusivity  of  iodide  ions  to  sulfate  ions  of  21.  Since 
the  chloride  diffusion  coefficient  for  cement  paste  is 
essentially  identical  to  that  of  iodide  diffusion  [10]  the  ratio 
of  the  effective  chloride  diffusion  coefficient  (Dec)  to  sulfate 
diffusion  coefficient  (Dgs)  also  should  be  close  to  21. 

Unfortunately,  the  water  to  cement  ratios  (w/c)  of  the  cement 
paste  and  mortar  were  not  reported.  Since  no  water-reducers  were 
reported  to  be  added,  it  is  likely  that  the  w/c  ratios  were  above 

0.4. 

A major  benefit  of  the  radioactive  tracer  method  for  measuring 
the  depth  of  penetration  of  sulfate  into  cement-based  materials 
is  that  the  calcium  sulfate  added  as  a set-regulator  to  portland 
cements  is  not  radioactive  and  thus  will  not  be  detected. 

2.3  Determination  of  Sulfate  and  Chloride  Penetration  Profiles 
by  EDXA. 

Penetration  profiles  of  chloride  and  sulfate  ions  into  concrete 
were  studied  by  Feldman,  et  al.[ll]  by  using  energy  dispersive  X- 
ray  analysis  (EDXA) . Concretes  were  made  from  two  portland 
cements,  Canadian  Types  10  and  50,  respectively  equivalent  to 
ASTM  Type  I and  Type  V portland  cements,  and  slag  and  silica  fume 
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mineral  admixtures.  Data  were  reported  for  one  year  of  diffusion 
studies.  Depth  of  penetration  was  taken  to  be  the  depth  at  which 
the  concentration  of  diffusing  sulfate  or  chloride  ions  became 
equal  to  the  intrinsic  concentrations  of  the  concrete,  based  on 
chemical  analysis. 

Four  w/c  ratios,  0.35,  0.42,  0.5,  and  0.60  were  used,  and  denoted 
as  mixes  Ml,  M2,  M3, and  M4,  respectively.  Three  concrete 
systems,  SI,  S2,  and  S5,  were  made  at  each  w/c  ratio  as  follows: 
SI  - Type  10  Portland  cement;  S2  - 90%  Type  50  Portland  cement  + 
10%  silica  fume;  S5  - 75%  blast  furnace  slag  + 22%  Type  50 
Portland  cement  + 3%  silica  fume. 

Data  were  available  in  tabular  and  graphic  form  for  the  depth  of 
ingress  of  sulfate  and  chloride  ions  at  specific  times.  The 
conventual  method  of  analyzing  such  data  is  to  use  the  classical 
error  function  solution  [ 12 ] : 


C- 

- Ci, 


= ezfc 


X 

(4D,t) 


(1) 


Where  is  the  surface  concentration  of  either  chloride  or 
sulfate  ion, 

Cb  is  the  background  (or  initial)  concentration  of  either 
chloride  or  sulfate 

C is  the  concentration  eit  a specific  depth  and  time, 

X is  the  depth  where  the  concentration  of  diffusing  ions 
is  0 , i . e . , C = Cb , 

Dg  is  the  effective  diffusion  coefficient,  and 
t is  time. 

Another  approach  to  calculation  De  without  the  use  of  the  error 
function  has  been  provided  by  Bazant  [13].  The  advantage  of 
this  method  is  that  the  calculation  of  De  can  be  easily 
computerized.  He  derived  the  following  parabolic  equation: 

^ ^ = [1  - (1/3^2  ■ x/2D^t)^'^]^  <2) 


Solutions  for  eq.  (2)  are  consistent  with  those  for  eq.  (1)  over 
the  range  of  the  dimensionless  variable  x/2(Det)^^^  from  0 to 
differing  in  this  range  by  no  more  than  0.033. 
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Setting  the  concentration  C equal  to  the  background  concentration 
Cb  in  eg.  (2)  and  rearranging: 

Z?^  = xVl2t  (3) 

The  method  used  here  is  consistent  with  the  method  of 
dimensionless  groups  given  by  Pommersheim  and  Clifton  [14].  Eq. 
(3)  predicts  a dimensionless  time  ratio  of  tDg/x^  = 1/12.  Since 
this  ratio  is  much  less  than  unity,  the  conditions  are  far  from 
steady  state.  It  is  thus  appropriate  to  use  eq.  (3)  to  obtain 
approximations  of  the  effective  diffusivities  of  ions  in  concrete 
or  cement  paste. 

Using  eq.  (3)  the  diffusion  coefficients  were  calculated  for  the 
three  batches  and  four  mixes,  and  the  results  are  given  in  Table 
1.  The  results  are  consistent  with  those  reported  by  Spinks,  et 
al  [7],  since  in  all  cases  the  individual  chloride  ion  diffusion 
coefficients  were  higher  than  those  for  sulfate.  The  ratio 
Dec/Des  10  or  higher  except  for  the  w/c  of  0.6  with  the  S5 

concrete  system.  It  should  be  noted  that  the  cement  mixture  in 
S5  largely  ccmsisted  of  blast  furnace  slag,  while  systems  SI  and 
S2  were  mostly  portland  cement  (S2  contained  10%  silica  fume) . 


3.  MODELING  OF  SULFATE  DIFFUSION 

Two  approaches  for  modeling  sulfate  diffusion  are  described:  a 
finite  difference  model;  and  a "blind  ant  model."  These 
approaches  are  both  being  used  in  modeling  studies  at  the 
National  Institute  of  Standards  and  Technology  on  transport 
properties . 

3.1  Finite  Difference  Model  for  Sulfate  Diffusion/Reaction 

This  approach  used  a one-dimensional  model  for  the  ingress  of 
sulfate  ions  from  a constant  concentration  source  into  a slab  of 
hardened  cement  paste.  The  cement  paste  is  initially  free  of 
unbound  sulfate  ions  (C=0) . The  sulfate  ions  diffuse  into  the 
cement  and  are  consumed  in  an  irreversible  reaction  with  calcium 
hydroxide  acording  to: 

S0«‘  + Ca(0H)j  ->  CaSO^  * 20H-  (4) 
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Table  1.  Ratio  of  Chloride  Diffusivity  to  Sulfate  Diffusivity* 


Concrete  System:  SI 


w/c: 

0.35 

0.42 

0.50 

0.60 

(mVs) 

4.8x10'^^ 

6.6x10’^^ 

1.9x10’^2 

3.6x10"^^ 

D,s  (®Vs) 

4.9x10'^^ 

4.0x10'^^ 

1.3X10'^^ 

0.98X10'^^ 

Ratio  Dec/D^s 

10 

16 

15 

27 

Concrete  System:  S2 

w/c: 

0.35 

0.42 

0.50 

0.60 

Dec  (mVs) 

2.9x10'^^ 

3.4x10"^^ 

4.9x10'^^ 

1.6x10'^^ 

D,s  (mVs) 

1.9x10*^^ 

2.2X10’^^ 

1.7x10'^^ 

1.7xl0‘“ 

Ratio  D^/D,s 

15 

14 

29 

39 

Concrete  System:  S5 

w/c: 

0.35 

0.42 

0.50 

0.60 

Dec  (mVs) 

4.2x10'^^ 

1.9X10*^^ 

4.4x10'^^ 

1.1x10'^^ 

D,s  (mVs) 

2.7x10"^^ 

1.7X10'^^ 

4.5x10'^^ 

2.1x10'^^ 

Ratio  Dec/Des 

16 

11 

10 

5.2 

*Based  on  penetration  data  from  reference  no. 11. 


According  to  Crank  [12],  the  standard  differential  equation  for 
this  problem  is: 


dC 

dt 


= D 


d^C 
^ dx^ 


- kC 


(5) 
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where  C is  the  concentration  of  sulfate  ions  at  time  t and 
position  X,  Di  is  the  intrinsic  diffusion  coefficient  (mVs)  and  k 
is  the  reaction  rate  constant  (s'^) . The  development  assumes 
that  sulfate  is  consumed  by  a first  order  chemical  reaction  and 
that  transport  occurs  only  in  the  x direction,  i.e.,  one- 
dimensional diffusion.  It  also  assumes  that  the  concentration  of 
sulfate  ions  controls  reaction  (4)  above  since  the  concentration 
of  calcium  hydroxide  does  not  appear  in  eq.  (5) . The  solution  to 
eq.  (5)  for  this  system  is  also  given  by  Crank  [12]: 


_C 

Co 


=-iexp(-x,  -^)erfc{ — — — -JI^) 
2 ^ N 2yQ5Jt 


+~exp  {x^ 


..  A)erfc{-^^^yPu) 
N 2,/UTE 


(6) 


Eq.  (6)  could  be  plotted  to  provide  concentration  vs.  position 
for  any  given  time  t.  The  difficulty  arises  in  that  an  infinite 
supply  of  calcium  hydroxide  does  not  exist  in  the  cement  paste  so 
that  at  each  position  in  the  sample  of  interest,  reaction  (4) 
will  terminate  at  some  specific  time. 

Since  this  constraint  cannot  be  conveniently  handled 
analytically,  a finite  difference  method  is  applied  to  this 
problem.  The  spatial  and  temporal  derivatives  in  eq. (6)  are 
approximated  by: 


/ yC  \ _ / ^m*l  \ 

0x2^"*  ‘ (Ax)  2 

t dC\  _ 

At 


(7) 


where  subscripts  indicate  spatial  locations  and  superscripts 
indicate  points  in  time  respectively.  Ax  and  At  indicate  the 
step  sizes  in  space  and  time  for  the  finite  difference  method. 
The  equations  in  (7)  can  be  substituted  into  eq.  (5)  and  giving 
an  iterative  solution  for  the  new  concentration  at  location  m: 


c:  = 


(Ax) 


(8) 


where  is  a flag  (with  values  of  0 or  1)  indicating  if  reaction 
can  occur  at  position  m.  is  set  equal  to  0 when  all  of  the 
calcium  hydroxide  is  consumed  at  a given  depth  m.  At  the  far  end 
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f 
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(zero  flux)  boundary  condition  is  used  so  that  eq. (8)  becomes 
[14]: 


C;  = + At[— ^(2C^i-2Cj  - kCJ(J 

(Ax^) 


(9) 


Eqs.  (8)  and  (9)  can  be  solved  using  the  computer  to  give  sulfate 
concentration  as  a function  of  time  and  position.  An  initial 
study  was  conducted  with  the  following  parameters: 


1.0  moles/ liter 
0.01  liters 
0.001  meters 
1*10'^^  mVs 
0,1*10“%  1*10’®  s“^ 

0.22 

1*10^  seconds 

50000  or  5*10®  seconds 


Surface  concentration- 
volume  of  a layer- 
Thickness  of  a layer  (ax)- 
Dif fusion  coefficient  (Di)- 
Reaction  constant  (k) - 
CH  volume  fraction- 
Time  step  (At)- 
Total  number  of  steps- 


The  calcium  hydroxide  volume  fraction  of  0.22  would  be  typical 
for  a fully  hydrated  cement  paste  with  a w/c  ratio  in  the  range 
of  0.4  to  0.5  [15]. 

The  results  are  shown  in  Fig.  1 for  times  of  1x10®,  5x10®,  and 
5x10®  seconds,  corresponding  to  about  3.2,  15.9,  and  159  years 
respectively.  The  effects  of  reaction  on  impeding  the 
penetration  of  sulfate  ions  into  the  cement  paste  are  evident 
from  the  plots  in  Fig.  1.  The  value  of  the  reaction  constant  k 
had  a significant  effect  on  the  progress  of  sulfate  penetration. 
Thus,  this  parameter  and  the  diffusion  coefficient  must  be 
known  to  predict  the  penetration  of  sulfate  ions  into  a cement 
paste  or  concrete  specimen. 


3.2  Blind  Ant  Model  for  Sulfate  Diffusion  in  Cement  Paste 

An  alternate  approach  for  simulating  diffusion/reaction  processes 
in  a random  porous  media  is  the  utilization  of  random  walkers  as 
described  in  detail  in  [16].  This  approach  was  first  applied  to 
diffusion  problems  in  the  early  1950s  [17].  Here,  walkers 
(commonly  known  as  blind  ants)  are  released  at  random  locations 
in  the  pore  space  of  a porous  media.  The  walkers  undergo  random 
Brownian  motion-type  diffusion  with  the  restriction  that  they 
cannot  enter  the  solid  phase  of  the  porous  media.  When  a step 
would  land  a walker  in  the  solid  phase,  the  step  is  not  allowed, 
but  the  timer  still  advances.  By  monitoring  the  average  distance 
a walker  travels  over  time,  a diffusion  coefficient  for  the 
porous  media  can  be  determined  [ 18 ] . 
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This  basic  approach  was  modified  to  study  the  diffusion/ reaction 
of  sulfate  ions  in  cement  paste.  First,  instead  of  releasing 
random  walkers  throughout  the  pore  space,  a constant 
concentration  of  walkers  was  maintained  at  one  surface  of  the 
porous  media  to  simulate  the  experimental  boundary  condition. 
Second,  the  walkers  were  allowed  to  diffuse  freely  in  both  the 
capillary  pore  space  and  the  nanoporous  calcium  silicate  hydrate 
(C-S-H)  gel  phase.  Finally,  as  the  walkers  are  diffusing,  they 
may  be  consumed  in  a reaction  at  the  surface  of  a calcium 
hydroxide  crystal  according  to  the  stoichiometry  given  in 
equation  (4) . 

This  approach  has  been  used  previously  to  study  the 
diffusion/binding  of  chloride  ions  in  hardened  cement  paste  [18]. 
To  do  this,  the  cement  paste  microstructure  is  represented  in  the 
computer  as  a two-dimensional  array  of  elements  (called  pixels) . 
Each  element  is  uniquely  identified  as  one  of  the  phases  of 
cement  paste;  in  its  simplest  form  the  phases  to  be  considered 
are  capillary  porosity,  unreacted  cement,  C-S-H  gel,  and  calcium 
hydroxide.  This  two-dimensional  representation  could  be  obtained 
from  a scanning  electron  microscope  image  of  a real  hydrated 
sample  or  from  a digital-image-based  computer  simulation  [19]. 
Fig.  2 shows  a 250x1000  digital  image  created  using  the  NIST 
cement  hydration  microstructural  model  [18,19].  Each  element  is 
about  1 in  size  so  that  the  sample  shown  in  Fig.  2 would  be 
0.25  by  1.0  mm.  White  indicates  calcium  hydroxide,  light  grey  is 
unreacted  cement,  dark  grey" is  C-S-H  gel,  and  black  is  porosity. 
The  original  system,  having  a w/c  ratio  of  0.5,  was  created  from 
a data-base  of  real  cement  particle  shapes  obtained  by  a 
combination  of  scanning  electron  microscopy  and  image  analysis 
[20].  In  the  hydrated  image  shown  in  Fig.  2,  a degree  of 
hydration  of  0.87  was  achieved  using  the  microstructural  model. 

This  image  can  be  used  as  input  into  a computer  simulation  for 
the  diffusion/reaction  of  sulfate  ions  in  cement  paste.  As 
outlined  previously,  the  sulfate  ions  are  represented  by  random 
walkers  which  can  diffuse  from  pixel  to  pixel  within  the  porosity 
and  C-S-H  gel  phases  and  react  with  the  solid  calcium  hydroxide, 
forming  calcium  sulfate,  according  to  eg.  (4) . By  initially 
placing  and  subsequently  maintaining  a constant  number  of  random 
walkers  at  the  top  surface  of  the  system,  diffusion  under  a 
concentration  gradient  can  be  simulated.  The  reaction  rate  of 
calcium  hydroxide  with  the  sulfate  ions  can  be  controlled  by 
designating  a probability  of  reaction  when  a diffusing  sulfate 
species  "collides”  with  a solid  calcium  hydroxide  pixel.  By 
specifying  this  probability  to  be  zero,  the  no  reaction  case  can 
be  simulated  as  a control  system.  By  tabulating  how  many  sulfate 
species  are  consumed  by  each  pixel  of  calcium  hydroxide,  the 
local  reaction  can  terminated  when  all  of  the  solid  calcium 
hydroxide  at  a specific  pixel  has  been  depleted  according  to  the 
stoichiometry  in  the  reaction  of  eg.  (4) . As  in  the  finite 
difference  solution,  the  concentration  of  sulfate  species  as  a 


8 


function  of  depth  can  be  monitored  throughout  the  simulation  to 
produce  the  concentration  vs.  distance  curves  counterpart  to 
those  in  Fig.  1. 

Since  each  sulfate  walker  can  represent  a collection  of  diffusing 
ions,  the  simulation  was  executed  such  that  one  walker  per  pixel 
corresponded  to  a concentration  of  0.1  M.  Maintaining  2500 
walker  along  the  top  surface  of  the  microstructure  results  in  a 
surface  concentration  of  1 M.  Results  for  this  simulation  are 
shown  in  Fig.  3 for  three  different  reaction  probabilities  (k) 
after  200000  steps  have  been  teJcen.  The  form  of  the  results  is 
seen  to  closely  follow  that  of  those  shown  in  Fig.  1,  validating 
that  the  "blind  ant"  method  can  produce  results  equivalent  to  the 
finite-difference  approach.  The  advantage  of  the  blind  ant 
method  is  that  it  can  be  applied  to  a specific  heterogeneous 
cement  paste  microstructure  (real  or  simulated) , whereas  the 
finite  difference  approach  considers  the  microstructure  to  be 
homogeneous  or,  at  most,  to  consist  of  a series  of  homogeneous 
layers . 

Both  techniques  clearly  illustrate  that  reaction  can 
significantly  slow  the  ingress  of  sulfate  ions  into  a cement 
paste  specimen.  To  quantify  this  effect,  the  kinetics  of  the 
reactions  occurring  must  be  measured  either  directly  or 
indirectly.  As  Fig.  1 shows,  the  kinetics  could  be  inferred  by 
monitoring  the  advance  of  a "sulfate  front"  into  a specimen  over 
time.  The  methods  described  in  section  2 could  be  used  to 
monitor  the  movement  of  a sulfate  front  or,  if  not  present,  to 
monitor  the  maximum  depth  of  ingress  of  sulfate  ions  as  a 
function  of  time. 

The  simulation  models  described  in  this  paper  should  be  further 
developed  to  provide  a means  of  extending  experimental  studies  on 
the  diffusivity  of  both  chloride  and  sulfate  ions  in  concrete  to 
other  concrete  mixtures  and  exposure  conditions.  Sulfate 
diffusion  studies  are  especially  difficult  to  carry  out  and 
with  further  development  it  is  reasonable  to  expect  that 
modeling  should  provide  a means  for  estimating  sulfate  diffusion 
coefficients. 


4.  DISCUSSION 

4.1  Comparison  of  Diffusivity  of  Chloride  and  Sulfate  Ions 

The  diffusion  coefficient  of  chloride  ions  and  sulfate  ions  at 
infinite  dilution  in  bulk  water  (free  stream  values)  are 
2.03x10^  and  1.07x10’*  ®^/s,  respectively  [21].  These  values  are 
some  three  to  five  orders  of  magnitude  larger  than  the 
^iif fusivities  in  cement-based  materials  because: 

(1)  the  effective  area  of  diffusion  is  lowered  by  the 
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presence  of  the  solid  phase. 

(2)  the  connectivity  of  the  pores  may  not  be  continuous 

(3)  the  diffusion  path  is  tortuous  which  increases  the  path 
length 

(4)  a given  diffusion  path  has  constricted  regions  which 
give  a greater  resistance  to  the  diffusive  flow  of 
molecules  or  ions. 

The  effects  of  these  constraining  factors  are  mathematically 
expressed  in  the  relationship  [22]: 


(11) 

The  porosity  c (volume  fraction  pores)  and  constriction  factor  o 
are  both  less  than  unity,  while  the  tortuosity  factor  x is 
greater  than  unity.  These  factors  act  to  lower  the  intrinsic 
diffusion  coefficient,  D^f^,  of  species  A relative  to  its  free 
stream  value  Because  the  geometrical  factors,  o and  x are 

complex  functions  of  pore  structure,  it  is  not,  at  present, 
feasible  to  use  eg.  (11)  for  reliable  calculations  of  DeA* 

Another  factor  which  can  reduce  the  intrinsic  diffusion 
coefficient  which  is  not  often  considered,  is  the  effect  of 
concentration.  Increasing  the  concentration  of  diffusing  ions 
will  tend  to  decrease  [12]^ 

The  diffusion  coefficient  of  chloride  and  sulfate  ions  diffusing 
through  the  same  cement-based  materials  will  both  be  equally 
effected  by  the  above  physical  factors.  Thus,  the  ratios  of  the 
intrinsic  diffusion  coefficients,  Dic/Dis,  should  be  approximately 
2.  However,  this  analysis  does  not  consider  the  effects  of 
binding  by  reaction  and  absorption  processes.  The  ratios  for 
effective  diffusion  coefficients,  Dgc/Des,  given  in  section  2, 
indicate  that  sulfate  diffusion  is  significantly  more  strongly 
affected  by  binding  than  chloride  ions.  Chloride  ions  are  known 
to  react  with  the  hydration  products  of  tricalcium  aluminate  and 
the  ferrite  phase  (4Ca0-Al203-Fe0)  in  hydrated  cement  to  form 
Friedel’s  salt  (3CaO-Al2O3-CaCl2-10H2O)  [23];  whereas,  sulfate 
ions  penetrating  into  concrete  are  likely  to  be  bound  by  reaction 
with  calcium  hydroxide,  tricalcium  aluminate  hydration  products 
such  as  monosulfate  (3Ca0-Al203-CaS04-12H20) , and  unreacted 
tricalcium  aluminate  [24].  Because  there  is  more  calcium 
hydroxide  present  in  hydrated  portland  cement  than  aluminates, 
it  is  reasonable  to  assume  that  the  effective  diffusion 
coefficient  of  sulfate  ions  would  be  less  than  that  of  chloride 
ions.  In  making  this  analysis,  it  is  recognized  that  the  binding 
capacity  of  a concrete  for  chloride  and  sulfate  ions  may  be 
dependent  on  the  its  mix  proportions  and  constituents,  especially 
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if  mineral  admixtures  are  added. 


4.2  Predicting  the  Diffusion  Coefficient  for  Sulfate  Ions 

An  approach  for  estimating  the  diffusion  coefficient  of  sulfate 
ions  in  concrete  is  demonstrated  in  the  following.  The  diffusion 
coefficient  of  chlorides  as  a function  of  the  w/c  ratio  of 
concrete  can  be  calculated  using  equations  such  as  [25]: 

D^jim^/s)  = 3.17 xIO-^M2460( w/c  - 0.30)2  + 15]  (12) 


Based  on  the  analysis  in  Table  1 and  the  study  by  Spiles,  et  al. 
[9],  a conservative  value  of  0.10  is  used  for  the  ratio  of 
Des/Dec*  Introducing  this  ratio  into  eg.  (12)  , the  equation  for 
sulfate  diffusion  coefficient  becomes: 


Dso4{ni^/s)  = 3.17x10-^^2460  (w/c  - 0.30)  2 + 15]  (13) 


For  exeunple  if  the  w/c  ratio  is  0.4,  then  the  predicted  sulfate 
diffusion  coefficient  is  1.3x10"^^  m^/s.  Note  that  eqs.  (12)  and 
(13)  are  not  valid  if  w/c  < 0.3,  otherwise  they  would  predict 
that  the  effective  diffusion  coefficient  would  increase  as  w/c 
decreases  below  0.30. 

The  relationship  between  and  w/c  in  eq.  (12)  was  used  for 
demonstration  purposes.  Other  equations  have  been  reported  [26] 
for  calculating  chloride  ion  diffusion  coefficients  as  a function 
of  w/c  ratios  that  are  equally  credible  as  is  eq.  (12) . They  can 
be  modified,  as  described  above,  to  estimate  diffusion 
coefficients  for  sulfate  ions  intruding  concrete. 

4.3  Effect  of  Sulfate  Diffusion  Coefficient  Values  on  Service 
Life  Predictions 

Two  different  types  of  models  have  been  developed  for  making 
predictions  of  service  life  when  the  controlling  degradation 
process  is  sulfate  attack:  1)  a shrinking  core  model,  and  a 2) 
reaction  controlled  model  (or  mechanistic  model) . 

In  the  shrinking  core  model  [27],  which  is  diffusion  controlled, 
sulfate  ions  move  through  degraded  concrete  to  the  interface  of 
unreacted  concrete,  and  then  rapidly  reacts  with  the  hydration 
products  of  tricalcium  aluminate  to  form  expansive  products  such 
as  ettringite.  Mass  transport  equations  are  used,  assuming  a 
quasi-steady  state,  to  predict  the  movement  of  sulfate  within  the 
concrete.  The  flux  of  sulfate  ions,  N,  is  given  by: 
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N = D^{Cjx) 


(14) 


where  Co  is  the  concentration  of  sulfate  in  the  bulk  solution, 
is  the  effective  diffusion  coefficient,  and  x is  the  depth  of 
attack . 

The  rate  of  attack,  dx/dt,  is  given  by  the  rate  of  mass  transport 
divided  by  the  concentration  of  solid  tricalcium  aluminate,  C^: 

dx/dt  = -N/C^  = (15) 

Integrating  eg.  (15)  gives: 


(16) 


If  X is  the  depth  of  sulfate  attack  which  results  in  unacceptable 
performance  of  a concrete  element  or  structure,  then  according  to 
eq.  (16)  , a ten-fold  increase  in  will  decrease  the  service  life 
by  10 . 

In  the  reaction  controlled  model  [3],  it  is  assumed  that: 

(1)  sulfate  ions  from  the  environment  penetrate  the 
concrete,  by  diffusion, 

(2)  sulfate  ions  react  expansively  with  aluminates 
in  the  concrete,  and 

(3)  cracking  and  deleunination  of  concrete  surfaces  result 
from  the  expansive  reactions. 

Cracking  and  delamination  of  the  concrete  surface  exposes  a new 
surface  to  a concentration  of  sulfate  ions  near  that  of  the 
groundwater  sulfate  concentration,  rather  than  the  smaller 
concentration  resulting  from  diffusion.  The  model  predicts  that 
rate  of  sulfate  attack  will  be  largely  controlled  by  the 
concentration  of  sulfate  ions  and  aluminates,  diffusion  and 
reaction  rates,  and  the  fracture  energy  of  concrete.  The 
degradation  rate  (R)  is  linear  in  time  (m/s)  and  is  given  by  [3}: 

R “ Xapaii/Tapaii  = ( EB^CgCoeDi ) / ( ai|F  (1”V)  (17) 

where  Xspan  is  the  thickness  of  a spalled  layer, 

Tspaii  is  the  time  for  a layer  to  spall, 

E is  Young's  modulus, 

B is  the  linear  strain  caused  by  one  mole  of  sulfate. 
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reacted  in  1 of  concrete, 

Cs  is  the  sulfate  concentration  in  bulk  solution, 

Coe  is  the  concentration  of  reacted  sulfate  as  ettringite, 

Di  is  the  intrinsic  diffusion  coefficient  of  sulfate  ions, 
tto  is  a roughness  factor  for  fracture  path  (assumed  to  be 
1.0)  , 

i|r  is  the  fracture  surface  energy  of  concrete,  and 
V is  Poission's  ratio. 

Note  that  both  the  shrinking  core  model  and  the  reaction  controlled 
model  predict  that  a 10-fold  increase  in  the  diffusion  coefficient 
of  sulfate  will  result  in  a 10-fold  decrease  in  service  life. 
Therefore,  representative  estimates  of  sulfate  diffusion 
coefficients  are  needed  to  madce  reasonable  service  life 
predictions . 


5.  SUMMARY  AND  CONCLUSIONS 

The  sulfate  diffusion  coefficient  is  a variable  in  theoretical 
models  used  for  predicting  the  service  life  of  concrete.  Models 
based  on  diffusion  control  and  reaction  control  predict  that  a 10- 
fold  increase  in  De  will  decrease  the  predicted  service  life  by  a 
factor  of  10.“  Therefore,  reliable  estimates  of  sulfate  diffusion 
coefficients  are  needed  to  make  reasonable  service  life 
predictions . 

Methods  for  determining  the  sulfate  diffusion  coefficient  of 
sulfates  in  cement-based  memorials  are  explored,  including 
experimental  methods  and  mathematical  modeling.  In  experimental 
studies  the  depth  of  penetration  of  sulfate  ions  has  been  monitored 
energy  dispersive  x-ray  analysis  (EDXA)  techniques.  A major 
benefit  of  the  radioactive  tracer  method  is  that  the  calcium 
sulfate  added  to  the  clinker  will  hot  be  radioactive  and  thus  not 
be  detected. 

Sulfate  diffusion  was  modeled  by  two  approaches:  a finite 
difference  model;  and  a "blind  ant  model.”  Both  modeling 
approaches  gave  similar  results.  The  blind  ant  method  has  the 
advantage  that  it  can  be  applied  to  a specific  heterogeneous 
microstructure,  while  the  finite  difference  method  is  one- 
dimensional and  must  consider  the  microstructure  to  be  homogeneous 
or  at  most  to  consist  of  a series  of  homogeneous  layers.  Both 
models  show  that  the  ingress  of  sulfates  into  concrete  is  highly 
dependent  on  both  its  diffusion  coefficient  and  reaction  constant. 
Further  development  of  the  models  is  necessary  before  they  can  be 
used  to  obtain  reliable  predictions  of  diffusion  coefficients. 

/ui  approach  for  estimating  the  diffusion  coefficient  of  sulfate 
ions  has  been  developed  which  is  based  on  comparing  sulfate 
diffusion  coefficients  with  chloride  diffusion  coefficients.  For 
concretes  in  which  portland  cement  is  the  predominant  cementitious 
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material,  a Dec/Des  ratio  of  10  was  proposed.  This  value  can  be 
combined  with  equations  which  predict  chloride  ion  diffusivity  as 
a function  of  w/c  ratio,  to  estimate  the  sulfate  diffusion 
coefficient. 
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8.  NOMENCLATURE 


B 


C 

Cb 

Co 

Coe 

C3 

C3 

De 

Dec 

DeS 

Di 

DiA 

E 

k 

m 

N 

K 

t 

^apall 

V 

w/c 

X 

^spall 


- linear  strain  caused  by  one  mole  of  sulfate 
reacted  in  1 m^  of  concrete 

- concentration  at  specific  depth  and  time 

- background  or  initial  concentration 

- concentration  at  position  m 

- concentration  of  sulfate  ions  in  bulk  solution 

- concentration  of  reacted  sulfate  as  ettringite, 

- surface  concentation 

- sulfate  concentration  in  bulk  solution 

- concentration  point  in  time 

- effective  diffusion  coefficient 

- effective  diffusion  coefficient  for  chloride  ions 

- effective  diffusion  coefficient  for  sulfate  ions 

- intrinsic  diffusion  coefficient 

- intrinsic  diffusion  coefficient  of  species  A 

- Young's  modulus 

- reaction  rate  constant 

- position 

- flux  of  sulfate  ions 

- flag 

- time 

- time  for  a layer  to  spall 

- Poission's  ratio 

- water-to-cement  ratio 

- depth 

- thickness  of  a spalled  layer 


At  - step  size  in  time 

AX  - step  size  in  space 

€ - porosity 

o - constriction  factor 

«o  ” roughness  factor  for  fracture  path 

X - tortuosity  factor 

♦ - fracture  surface  energy  of  concrete 
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CONCENTRATION  (M) 
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Fioure  1.  Concentration  of  diffusing  sulfate  in  cement  paste 

predicted  by  finite  difference  model.  Exposure  times  to 
1 N solution  solution  are  identified  by  A,  B,  and  C, 
having  respective  values  of  3.2,  15.9,  and  159  years. 

K is  the  reaction  constant. 


18 


i"--]  “'^^^4*'®  Si«i.8tfHib  5o  mlanutSmofmo  ,t  tt-umi'i 


Figure  2.  A 250x1000  digital  image  created  using  the  NIST  cement 
hydration  microstructural  model  [18,19].  White 
indicates  calcium  hydroxide,  light  grey  is  unreacted 
cement,  dark  grey  is  C-S-H  gel,  and  black  is  porosity. 
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Figure  3.  Concentration  of  diffusing  sulfate  in  cement  paste 
predicted  by  blind-ant  model.  P (React)  is  the 
probability  of  reaction. 
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